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Abstract —Wideband spectrum sensing is a significant challenge 
in cognitive radios (CRs) due to requiring very high-speed analog- 
to-digital converters (ADCs), operating at or above the Nyquist 
rate. Here, we propose a very low-complexity zero-block detection 
scheme that can detect a large fraction of spectrum holes from the 
sub-Nyquist samples, even when the undersampling ratio is very 
small. The scheme is based on a block sparse sensing matrix, 
which is implemented through the design of a novel analog-to- 
information converter (AIC). The proposed scheme identifies some 
measurements as being zero and then verifies the sub-channels 
associated with them as being vacant. Analytical and simulation 
results are presented that demonstrate the effectiveness of the 
proposed method in reliable detection of spectrum holes with 
complexity much lower than existing schemes. This work also 
introduces a new paradigm in compressed sensing where one is 
interested in reliable detection of (some of the) zero blocks rather 
than the recovery of the whole block sparse signal. 

Index Terms - Wideband spectrum sensing, cognitive radios, 
compressed sensing, block sparse signals. 

I. Introduction 

In a wireless communication environment, many of the 
primary users (PUs) do not use their licensed frequency bands 
at all times. The surveys show that the maximum frequency 
utilization of the allocated spectrum can be less than 10% 
111. To increase the frequency utilization in such environments, 
secondary users (SUs) equipped with cognitive radios (CRs) (2| 
can be deployed to use the available spectrum (spectrum holes) 
without interfering with the PUs. Spectmm sensing, as the task 
of finding the spectmm holes, is thus one of the important 
functions performed by CRs. 

As an integral part of spectrum sensing, an analog to digital 
converter (ADC) is often used to sample the signal at or above 
the Nyquist rate. In wideband scenarios, this will require very 
high-speed ADCs, which are challenging to implement. 

One of the techniques of wideband spectmm sensing (WSS), 
which could relax the high-speed sampling requirement, is 
a channel-by-channel scanning approach. A tunable bandpass 
filter (BPF) scans one channel at a time. A narrowband spec¬ 
tmm sensing technique is then applied to sense if the scanned 
channel is vacant or not. This approach is, in general, slow 
and inflexible Another approach, called “filter-bank” g), 
is based on using a parallel set of narrowband filters, each 
covering one channel. In this approach, however, an enormous 
number of RF components is required due to the parallel 
structure of the filter bank. 


Since the spectmm is sparsely occupied, the information 
rate of the analog signal is much less than the Nyquist rate. 
Therefore, an analog to information converter (AIC) can be 
used with sub-Nyquist sampling rates, based on the compressed 
sensing (CS) theory 0, 

Compressed sensing is a signal processing technique that can 
recover an unknown signal with high probability from a small 
set of linear projections, called measurements , if the signal is 
sparse in some domain 0. Linear projections of the signal and 
its sparse representation are described by two matrices called 
measurement matrix and sensing matrix, respectively. 

CS has been recently applied to the wideband spectrum 
sensing to alleviate the need for high sampling rates ©-HD- 
The spectrum is then reconstructed by a sparse signal recovery 
algorithm, see, e.g., To reduce the sampling rate, in 

(To). (TT), the block structure of the spectrum is also used in the 
reconstruction process[[] In all the above literature [81—|Tl|, the 
spectrum holes are detected from the reconstructed spectrum 
by applying an energy detection algorithm. Although, in each 
CR, only finding one vacant sub-channel may be enough, in 
cooperative CR networks, finding more than one spectrum hole 
is of interest. Therefore, in this paper, our goal is to detect 
as many spectrum holes as possible with a given number 
of measurements. Neither the recovery of the signal, nor the 
detection of all spectrum holes is of direct interest. We thus 
introduce a new CS paradigm in which, one is interested in 
reliable detection of some of the zero blocks, even when the 
number of measurements is relatively low. 

To the best of our knowledge, all wideband spectmm sensing 
methods in the literature that use CS for spectrum sensing 
(e.g., @, (TO), (TT)), are based on dense sensing matrices 
(matrices which have few, or no, zero entries). The recovery 
algorithms with such matrices are often based on linear or 
convex optimization with the complexity between 0(N 2 ) and 
0(N 3 ), where N is the dimension of the signal. In mobile 
CRs with limited power and computational resources, it is 
impractical to use CS with such costly computations for the 
reconstruction of the spectrum. In addition, the processing time 
required for signal reconstruction can impose a significant delay 
in the spectrum sensing. These motivate us to look for an ultra 
low-complexity scheme for spectrum sensing with the goal of 

'block sparse signals are the sparse signals that have nonzero entries 
occurring in clusters. 


detecting a large proportion of spectrum holes. 

A popular category of recovery algorithms for CS are those 
based on message-passing algorithms. In particular, the ap¬ 
proximate message-passing (AMP) algorithm of G3’ which 
works with dense sensing matrices, has attracted much attention 
due to its remarkable performance/complexity trade-off. In the 
context of block sparse signals, it is demonstrated in m 
that AMP with James-Steins shrinkage estimator (AMP-JS) 
can outperform the existing block sparse recovery algorithms. 
Another subcategory of iterative message-passing recovery al¬ 
gorithms are verification-based ones [|14|-fT6), which are based 
on sparse sensing matrices and are thus of lower complexity 
compared to the message-passing algorithms on dense matrices 
(graphs). These algorithms however are very sensitive to noise. 
In addition, their application to the recovery of block sparse 
signals face a number of fundamental challenges. 

In this paper, we apply an ultra low-complexity scheme to 
the wideband spectrum sensing, which is based on a block 
sparse sensing matrix. To this end, we design a measurement 
matrix as part of the AIC, such that the resulting sensing 
matrix has a block sparse structure. As will be seen later, 
the structure of the sensing matrix affects the performance of 
zero detection scheme. We analyze the proposed zero-block 
detection scheme for regular sensing matrices (those with the 
same number of non-zero blocks in each row and column) and 
then generalize the results to the case of irregular matrices. We 
can optimize the structure of the sensing matrix to maximize 
the fraction of spectrum holes that can be detected for a 
given number of measurements. We consider both noiseless 
and noisy scenarios in this work, and analyze the performance 
of the proposed algorithm in both cases. Our analytical and 
simulation results demonstrate the effectiveness of the proposed 
algorithm in reliable detection of spectrum holes (even when the 
number of measurements is relatively small) with complexity 
that is negligible compared to the existing block sparse recovery 
algorithms such as AMP-JS. 


II. Notations, Definitions and Ensembles of 
Sensing Matrices (Graphs) 

In this paper, we use capital bold letters to refer to matrices 
and lowercase bold letters for vectors. We also use C and R to 
refer to the fields of complex and real numbers, respectively. 
For the real and imaginary parts of a complex argument, we use 
5ft(.) and Q(.), respectively. In addition, (.) T , (.)*, (.) _1 and (,) f 
denote the transpose, conjugate, matrix inversion and conjugate 
transpose, respectively. 

The zero block detection scheme presented in this work can 
be described based on a bipartite graph representation of CS 
operation. Let Q(V u M, S VtM ) denote a bipartite graph where 
V U M is the node set and £ V]M is the edge set, so that, every 
edge in S v M connects a node in V to a node in M. We refer 
to the sets V and M as variable nodes (VNs) and measurement 
nodes (MNs), respectively. In a bipartite graph, £s 1 ,s 2 C £v,m 
denotes the set of the edges that connect the set Si C V to 
the set S 2 C M. In addition, in a bipartite graph, let V(m) and 
M(v) denote the set of VNs neighboring to the measurement 


node m and the set of MNs neighboring to the variable node 
v, respectively. 

Let V z C V, V„ 8 C V, M z C M and C M refer to the 
sets of zero VNs, non-zero VNs, zero MNs and non-zero MNs, 
respectively. Clearly, V„ 8 (M „ z ) is the complement of the set 
V z (M z ) with respect to V (M). 

In a bipartite graph, a node in V {M) has degree i if it is 
connected to i nodes in M (V). We use d(v) to denote the 
degree of v. Let the polynomials A(x) = J2, Vaf and p(x) = 
J2, f ix’, respectively, represent the degree distributions of VNs 
and MNs, where A, and />, denote the fraction of degree-? nodes 
in VNs and MNs, respectively. Clearly, A(l) = p( 1) = 1 and 
A|V| = p\M\ = \S v ,m\, where A = JA?A», P = Y. i i Pi and l-l 
denotes the cardinality of a set. 

If the degree of all nodes in M and V are d M and d v , 
respectively, the graph is called bi-regular ( regular , in brief). 
Otherwise, it is called irregular. In a regular graph, p(x) = x dM 
and \(x) = x dv . 

III. System Model and Problem Statement 

Suppose a cognitive radio receives and down-converts a 
wideband analog signal, r(t), with the bandwidth of W. We 
assume that r(t) occupies L consecutive, non-overlapping spec¬ 
trum bands, each with bandwidth equal to B = W/L Hz, 
referred to as sub-channels. The boundaries of the sub-channels 
are denoted by f 0 < f y < • ■ • < f L . 

Suppose that the vector r £ R v is the discrete representation 
of the analog signal r(t), where A corresponds to the Nyquist 
rate. Let x = Fr be the spectrum of r with only K « N non¬ 
zero elements and F is the A x A unitary Fourier matrix. Since 
r(t) is a real signal, the spectrum x is conjugate symmetric. 
Therefore, we only focus on the positive frequency elements 
of x, denoted by + x. 

In the context of spectrum sensing, the spectrum is block 
sparse, i.e., at any given time, each sub-channel is occupied 
with probability a << 1 independent of the other sub-channels. 
We refer to a as the sparsity ratio. Let u (i) be the ?th segment 
of x, which corresponds to the elements in the range [/, ,. ), 

and / c w = (ft + /j_i)/2 be the center frequency of the ?th sub¬ 
channel. Therefore, we can write + x = [(u (1) ) T , ■ ■ • , (u (i) ) T ] T . 
We also assume that the magnitude of non-zero elements of x 
follows a distribution, /. 

In the noiseless scenario, the measurement vector y £ R M 
can be mathematically written as, y = d>r = 4>F _1 x = ©x, 
where © is the sensing matrix and d? £ R A/XJV is the mea¬ 
surement matrix. Since © = 4>F _1 and $ is a real matrix, 
each row of © has conjugate symmetry. Let + © denote the 
M x A/2 matrix including the last A/2 elements of each row 
of 0. In this paper, we assume that © is a random block sparse 
matrix in row with L blocks, in correspondence with the blocks 
of + x, and that only a few randomly selected blocks out of L 
blocks are non-zero and the rests are zeroQ In Section [TV] we 
will explain how a sensing matrix with this structure can be 
constructed. Assume that the magnitudes of all elements of a 

2 The sparse + © keeps the complexity low and improves the performance 
of the recovery algorithm, as becomes evident later. 
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Fig. 1: (a) The general structure of an AIC, (b) Sampling waveforms 
of the AIC for a block sparse sensing. 


non-zero block " follow a distribution g, where ' 0 : Jj be 
the 1th block of the mth row of + 0. The mth element of y can 
be written as 

j/ r „=25R^ + 0«.u (,) V m = 1,..., M, (1) 

where denotes the inner product. Now, the problem is to 
detect the maximum number of spectrum holes from a given 
number of measurements M. Note that here, we are interested 
in very small values of M, i.e., M « L. Such values of M 
can be much smaller than typical values that are of interest in 
conventional CS frameworks where M is small in comparison 
with N (not L). 

The relationships in |T]i for m = 1,..., M, can be represented 
by a bipartite graph, referred to as a sensing graph , with L VNs 
and M MNs. In this model, a VN is a sub-vector of B = N/(2L) 
complex elements. 

IV. Design of an AIC with Block Sparse Sensing 
Matrix 

In Fig. [TJa), the general structure of an AIC with M parallel 
branches of mixers and integrators (BMIs) is shown ED . In this 
section, the objective is to generate the sampling waveforms 
4>i{t) of Fig. |TJa) to create a block sparse sensing matrix. 

Suppose, © m and d>,„ are the mth row-vectors of © and 
d>, respectively. We have ©j„ = (F -1 ) 1 '^,, Vm. Since, F is a 
unitary matrix and <f> is a real matrix, we will have 0.J„ = 
F&f, Mm. In other words, ©,„ is the conjugate of the Fourier 
transform of <f\„ . Therefore, © will have block sparse rows, 
when the rows of the measurement matrix are block sparse in 
the frequency domain. 

Fig. |T]b) demonstrates how to generate the sampling wave¬ 
forms i = 1,..., M, in Fig. [TJa) to have a block sparse 

sensing matrix. In Fig. |T|b), the input signal, p c (t), is a noise 
generator with a rate equal to or higher than W/L. Then, we 
filter the input signal by a low-pass (LP) filter with the cut¬ 
off frequency of W/2L to result in a baseband signal with 
frequency range of [— W/2L, W/2L], Finally, using the carrier 
frequencies / c (1) , / c (2 \ • • ■ ,/ c (i) , the signal is modulated. Now, 
we can generate M sampling waveforms using a linear 

combination of the resulted waveforms vi(t),l = 1 ,,L, based 
on a properly designed mapping. Note that a mapping with 


sparse linear combinations is needed to create a sparse sensing 
matrix. 

The mapping in Fig. |T] b) can be constructed by the sensing 
graph described in Section[III] Therefore, we can make a corre¬ 
spondence between the VNs and waveforms v t (t), l = 1,... ,L, 
and between the MNs and sampling waveforms m = 

1 ,... ,M, as follows 

<M*) = 

lev(m) 

where |V(m)| = d(m). Let 4> m (f) and vff) denote the Fourier 
transforms of 4> m (t) and vi(t), respectively. Therefore, = 

Sigv(m) i- e -’ which consists of the samples of 
would be block sparse with d(m) non-zero blocks in +0 m . 

In Fig.jTjb), L multipliers are required. We can generate the 
sampling waveforms using only M multipliers, if we modulate 
the output signal of the LP filter by cos(27r/ c (i) f), m = 

L ■ • • : V7. 


V. Proposed Spectrum Hole Detection Scheme 


A. Noiseless Scenario: Proposed Scheme and Its Analysis 
A.l. Proposed Scheme 


The following lemma is the foundation of the proposed zero 
block detection in the absence of noise. 

Lemma 1: If y m = 0 in (JTJ, the sub-vectors of ' x cor¬ 
responding to the non-zero blocks of + © m , are zero, with 
probability one. 

Proof: The proof simply follows from the continuity of at 


least one the distributions / or g, described in Section III 


Based on Lemma 1, the proposed scheme detects a zero 
block in the sensing graph, if that block is connected to a zero 
measurement node. We refer to this simple detection scheme 
for zero blocks (vacant sub-channels) as zero measurement 
detection (ZMD). In order to quantify the performance of the 
proposed method, we define two measures: 


(i) Probability of zero block detection ( P Z d)'■ the ratio of the 
number of zero blocks detected correctly to the total number of 
zero blocks that exist in the spectrum. 

(ii) Probability of wrong zero block detection ( Pwzo ): the ratio 
of the number of occupied blocks detected falsely as vacant to 
the total number of the blocks detected as zero. 


Since P WZ d indicates the level of interference to the PUs, 
it is an important measure from the standpoint of PUs. Note 
that in the noiseless case, based on Lemma 1, Pwzn = 0 and 
the objective is to design a sensing graph in order to maximize 


A.2. Analysis 

In this section, we first consider bi-regular sensing graphs. In 
such graphs, the MN m £ M z , if and only if, all d M neighboring 
VNs belong to V,. Therefore, the probability that m is zero is 
calculated by 

P (m £ Mf) = (1 — a) dM . (2) 

Therefore, E(\M Z \) = M (1 — a) dM . 




















On the other hand, the zero VN, v, is recovered by the 
ZMD method, if and only if, at least one of the neigh¬ 
boring MNs is in M z . By definition, P Z d is the ratio of 
the expected number of detected zero variables to the total 
number of zero VNs. Therefore, it is clear that P ZD is the 
probability that the zero VN, v, is recovered by ZMD, i.e., 
Pzn = P (\£v, Mz | > l|t> £ V z ) = 1 — P | = 0|tt £ V,). 

For calculating this probability, let the edge e connect a VN 
in V* to a MN in M. We first calculate the probability p 0 that 
the edge e connects a VN in V z to a MN in M ,. Let /, and h e 
be the tail node and the head node of the edge e in M and V, 
respectively. Therefore, we can write 


Po 


P(t e £ M Z \K £ Vz) 

P(fe £ M z ,h e € Vz) _ P(t e £ M z ) 
P (ft. £ Vz) “ P(/le <= Vz) 
E(\£ v , Mz \)/\£ v , M \ = E(\£ V ' Mz \) 
E(\£ Vx .m\)/\£ v .m\ E(\£ Vz , m \) " 


d M M (1 - a) dM 
d v L(l — a) 


(1 — a ) dM ~ 1 , 


d M E(\M z \) 

d v |Vz| 

( 3 ) 


where the last equality comes from the fact that in bi-regular 
graphs, d v = Md M /L. Thus, the probability that the zero 
variable node v does not have any connection to M z , i.e., all 
d v neighboring MNs belong to M nz , is (1 — p 0 ) dv ■ Therefore, 


Pzn = 1 - P (|£„,Aiz I =0|u £ V z ) 

= 1 - (1 - p 0 ) dv = 1 - (1 - (1 - a) dM ~ 1 ) dv . (4) 

We also can simply generalize the analysis of bi-regular 
sensing graphs to the case of irregular ones with the constraint 
that the VNs and MNs follow the degree distributions \(x) and 
p{x), respectively. After following the same steps as the ones 
for regular graphs, we can obtain P ZD as 

p zd = i-J2Mi-PoY, ( 5 ) 

where p' 0 = JT ip, (1 - a) i_1 /Ei *P<- 


B. Noisy Scenario: Proposed Scheme and Its Analysis 
B.l. Proposed Scheme 

In practice, the measurements are corrupted by noise, i.e., 
y = ©x + n, where n is an M x 1 vector whose elements are 
independent and identically distributed (i.i.d.) zero-mean Gaus¬ 
sian random variables with variance ojj, J\f(0,cr*). Therefore, 
the mth measurement can be written as 



Noise has a destructive effect on the application of ZMD 
method. In other words, if the noise-free measurement is equal 
to zero, its noisy version is not equal to zero with probability 
one. This will, in effect, disable ZMD method as described 
for noiseless scenarios. Here, we propose a test to detect zero 
measurements and we also show that in the high SNR regime, 
the test is translated to a thresholding approach, and in this 
regime, we obtain the optimum threshold, analytically. Using 
the thresholding technique, the problem is reduced to one 
similar to the noiseless scenario. 


We consider an irregular sensing graph. For the analysis, we 
assume that the real and imaginary parts of non-zero elements 
of the sub-vectors u w both come from a Gaussian distribution 
Af(0, a-.), in an i.i.d. fashion. We consider d(m) hypotheses Hi 
for i = 0,..., d(m), where Hi represents the scenario that the 
measurement m is connected to i non-zero blocks and d(m) — i 
zero-blocks. Therefore, the a priori probability of Hi is 

P(Hi) = ^™^a i {l-a) dlrn) - i , i = 0,...,d(m). (7) 

Furthermore, suppose that only one of the neighboring blocks 
of the mth MN (e.g., the jth block) out of d(m) blocks is non¬ 
zero and d{m) — 1 neighboring blocks are zero (i.e.. Hypothesis 
Hi is true). In this case, we have 

y m = 2K ( + ©“ ) .u w) ) +n m 

D 

= 2^»( + ©W).3?(u“ ) ) - 3^( + ©^ ^ ).3^(u‘ J, ) + n m , (8) 

i= 1 

where ' ©;)), and ' are the zth element of the sub-vector 
+ ©^ 5 and the ith element of the sub-vector u (j> , respectively. 
The probability density function of y„, under this assumption 
(Hi) is 


P(ym\Hi) = AT (0,4|| + ©^ ) ||!cr“ + o 2 n ). 


(9) 


Based on the structure of the AIC presented in Section IV 
ll + ®m > l |2 = cte, Vj. Without loss of generality, we assume 
that || + ©)i ) ||2 = 1. Therefore, P(y m \Hi) = JV (0, 4a z + <r“). 
Similarly, the probability density function of y m under the 
hypothesis Hi is given by 


P(y m \Hi) = A f (0,4 io 2 s + ovi) , i = 0, ..., d(m). 


GO) 


In ZMD, we need to determine whether the noise-free version 
of y m , i.e., y T r l J, is zero or not. If y r r ‘ r f is zero, H a is true with 
probability one. Otherwise, Hi for i = 1,..., d(m) is true ( H 0 
is true), i.e., H 0 is false. To make a decision on whether ?/":/ is 
zero or not, we use the likelihood ratio test (LRT). The LR is 
defined by 


A(y m ) 


P(Vm\H 0 ) 

P(y™|«o)’ 


(ID 


and the test is performed as follows: 


if A(y m ) < c : accept Ho; 

if A (y m ) > c : reject H 0 (accept Ho), (12) 


where c is a properly selected threshold. 

We can calculate P(y m \H 0 ) in <H) b y 

P(y m \n 0 ) = E P(Vm\H z )P(Hi). 

r(rto) , =1 

Therefore, using ([7]), ([TO]), ( [TT] ) and (jT3j. we have 

( ,i( r J )a < (l-a)‘ l(m) ~W(0,4ta. a + ^) 


(13) 




1 - (1 - a) d 


( 14 ) 
















In the high SNR regime, where a 2 » a), we can neglect 
o' 2 ,, against a 2 in ( p~4| i. We can thus simplify the decision for 'H„ 
(i.e., A (y m ) < c) to 

\Vm\ < C, (15) 


where P (t e G G V„ z ) = Pfa- 

For the denominator of ( fl9| ), we have 

P (\£v,M%UM™ | 2 l) = 1 — P (\£v,McuMf | = 0) 

= 1 - (1 - P (t e G Ml U M Z)) dy ,(21) 


where c! is determined by the target P W zd as described in the 
next section. In other words, if \y m \ < c', we decide yZ? = 0 
and if \y m \ > c', we decide y"J ^ 0. 

B.2. Analysis 

Based on the likelihood ratio test, the measurement values 
are partitioned into two regions: 

n 0 = {y m : \y m \ < c}, and TZ t = {y m : \y m \ > c'}. (16) 

Based on the correct or the erroneous detection of zero/non¬ 
zero measurements, we can partition the set of measurements 
M into four different sets: Ml, MZ, Ml z and M™ z . The set 
M™ represents the set of non-zero measurements which are 
detected erroneously as zero, M z = {y m : y m G TZ a ,i-io}- Sim¬ 
ilarly, Ml = {y m : y m G TlcHo}, MZ z = {y m ■ y m $ Tlo,'H 0 } 
and M° nz = {y m ■ y m ^ TlcPo}- In the following, we first 
consider the case where the sensing graph is regular. 

The probability of false alarm (for zero measurements) P FA 
is defined as the probability that a measurement y m is detected 
as zero while at least one of its neighboring VNs belongs to V„ ;2 , 
i.e., P fa = f llf K P(y\'Ho)dy. For the signal model considered 
here, by using (|Y3 ]i, we have 


Pfa — 


E-i! ( d f )«•(! - a)*"- 1 f y€1Zo P(y| H,)dy 

1 - (1 - a) d 

E?“ ( d f )«*d - «) dM - , erf( V2(4 ;; §+ ^ ) ) 


1 - (1 - a) 


d(m) 


(17) 


where erf(x) := (2/^/tF) f* e _t ” dt. 

Probability of detection (for zero measurements), P D , is 
defined as the probability that a zero measurement is detected 
correctly as zero, i.e., P D = f u€Ko P(y\Ha)dy. For our signal 
model, P D is given by 


P.^rf^). (18) 

Clearly, E(\Ml\) = MP(H 0 )Pd = M( 1 - a) dM P D and 
E(\MZ\) = MP(H 0 )Pfa = M{ 1 - P(H 0 ))Pfa = M( 1 - (1 - 
a) dM )P F A- 

Now, we obtain P W zn and P Z n of the VNs based on the P D 
and Pfa of the MNs. Based on the definition, we have 


P (y g V„*||£„,.Mf luvi™| 2 l) 

P (\£v,M%UM™ | > 1|u G V„ z ) P (v G V„z) 


P (\£v,M<iUM™\ 2 l) 


(19) 


Since, with probability one, = 0 for v G V nz , we can 

write 


POSSUMS-1 2 l|ueV M ) =P(|f. fJM y| 2 1|U6V„) 

= l-P(|£„,/<ty| = 0|uGV„.) 

= i-(i-P(t e eMZ\KeV nz )) dv , (20) 


where 


P(te g MIUMZ) = 


E(\Ml\) + E(\MZ\) 
M 


= (1 — a) dM P D + (1 — (1 — a) dM )P FA - (22) 

Since, P(» G V„ z ) = a, from ( [20] ) and ( [2T] >, Pwzd is calcu¬ 
lated by 

n i - (i - p*a)* v 

1 - (1 - (1 - a) d MP D - (1 - (1 - a) d M)P FA ) dv 

(23) 

The probability of zero detection, P Z n, can be calculated by 

PzD = P (\£v,M%UM™ | 2 l|v G V.) 

= 1 — P (\£v,m%um™\ — 0|v G V.) 

= 1 — (1 — P (f e G Ml U MZ\h c G V z )) dv , (24) 

where 

P(t e G Ml U MZ, h B G V z ) 


P{t e G Ml U MZ\K G V z ) = 


P (K G V z ) 


E(\£v z ,M% I) + E(\£v z ,M™ I) 1 

E(\£v,m\) (1 — a) 

E(,\M c z \)d M +E(\Mf\)Y,i^{(,d M -i)( d M )a»(l- a ) d M-VPCH 0 ) 


Md M ( 1_Q! ) 


(25) 


= (1 - a) M P D + (1 - (1 - a) dM ~ 2 )PFA. 

Therefore, 

Pzd = 1 — ^1 — (1 — a) dM ~ 1 PD — (l — (1 — a) dM ~ 1 ) Pfa^ V . (26) 

When the measurements are noiseless, P D = 1 and P h A = 0. 
With these values for P D and P FA , Pwzd in ( |23[ i will be zero 
and P ZD in ( |26| ) is reduced to Q. With a large value for the 
threshold c', all the measurements are detected as zero, which 
results in P D = l and P r A = 1. We see from Equations (23 i 
and ( |26[ ) that, in this case, P WZD = a and P ZD = 1, which are 
expected. 

For given values of d M , d v , a s and cr„, the threshold d can 
be chosen to provide a target P wzD . We can then obtain P ZD 
using this threshold level. It is important to note that while 
P D and Pfa depend on the signal model, Pwzd and P ZD do 


not. Therefore, Equations (231 and (26 1 can be applied to other 
signal models. 

The above results can be generalized to irregular sensing 
graphs (matrices). The derivation of P ZD and P WZD is similar 
to (p~9|)-(|26|). For irregular sensing matrices, we have 


Pwzd — 


ail — JV Ai (1 — Pfa) 


1 — ^ (l ~ Ppa ~ (Pd — Pfa) Pi(l — a ) J ) 


and 


= 1 - A, (l - Pfa - (Pn - Pfa) 


E - 'O' 1 


Ej 3Pi 

In noiseless scenario, when P D = 1 and P FA = 0, from the 
above equations, we have P WZD = 0 and P ZD reduces to 0. 




















C. Complexity of the Proposed Algorithm 

In terms of the complexity, the proposed scheme has very low 
complexity as it only needs to identify the zero measurement 
nodes and verify the adjacent variable nodes as vacant sub¬ 
channels. This is negligible in comparison with any of the 
existing algorithms for the recovery of block sparse signals. In 
particular, one of the most powerful algorithms for the recovery 
of block sparse signals, in terms of performance/complexity 
trade-off, is AMP with lames-Stein’s estimator (AMP-JS) (13). 
In terms of complexity, the complexity of AMP-JS is O(MN), 
which is much higher than the minimal complexity of the 
proposed scheme. 

VI. Simulation Results 

In this section, we present some simulation results to demon¬ 
strate the effectiveness of the proposed method for detecting 
the spectrum holes in both noiseless and noisy cases. We 
consider an N dimensional signal with the positive frequency 
range of [0 : JV/2], containing L non-overlapping channels 
of equal bandwidth 21. We use an IV-point unitary Fourier 
matrix to map the input signal from the time to the frequency 
domain. To model the block sparse spectrum, each block is 
independently selected to be non-zero with probability a (from 
positive frequency elements). We assign the real and imaginary 
parts of the elements of non-zero blocks from a Gaussian 
distribution Af(0, af) in an i.i.d. fashion. In our simulations, 
we set a s = 1. It is important to note that the choice of these 
continuous distributions does not affect the performance of the 
ZMD method, in either noiseless or noisy scenarios. 

A. Noiseless Measurements 

Fig. 0 shows P Z d , obtained both analytically and by simu¬ 
lations, versus d M for different sparsity ratios. In this figure, 
P ZD of ZMD is depicted for the regular graphs with L = 1000 
and M = 500. Thus, d v = d M / 2. As seen from this figure, 
simulation and analytical results match perfectly. Furthermore, 
for each sparsity ratio, the optimum d M is depicted in terms 
of maximizing P Z n ■ We see that for smaller sparsity ratios, 
larger d M is required to maximize P Z d and with increasing 
the sparsity ratio, the optimum value of d M decreases. In 
particular, for larger sparsity ratios, the best degree distribution 
is (d v ,d M ) = (1, 2). 

Fig. [3] shows P Z D versus number of measurements, M, for 
different algorithms in a case with L = 500 and a = 0.25. Note 
that the x-axis is in log scale. In this figure, we see P Z n of the 
regular graphs with d v = 1 when d M is obtained by d v L/M 
for each M, as well as the regular graphs with d M = 4 when 
d v is obtained by d M M/L for each M. We have also plotted 
the curve corresponding to a graph, where each measurement 
node is connected to only one variable node and each variable 
node is connected to at most one measurement node. This case 
is denoted by “ltol” in Fig. [3] 

In Fig. [3] we have also plotted P ZD of AMP-JS, which is 
known as the-state-of-the-art in block sparse signal recovery, 
for two different block lengths 3 and 5. For AMP-JS, the 
sensing matrix is selected to be a dense matrix randomly 



Fig. 2: Pzd obtained by both analysis (dashed lines) and simulations 
(solid lines) versus dM for different sparsity ratios for bi-regular 
graphs with L = 1000 and M = 500. The stars are the maximum 
Pzd for each sparsity ratio. 



Fig. 3: Pzd versus M for a signal with a = 0.25 and L = 500, for 
different algorithms and sensing graphs. 

constructed with i.i.d. elements, where each element is a zero- 
mean Gaussian random variable with variance 1/M. In the 
AMP-JS, there is a thresholding function based on the t 2 - 
norm of each block. Thus, the zero and non-zero blocks can 
be detected easily. Since in the noiseless scenario, unlike the 
case for the proposed ZMD scheme, the P W zd of AMP-JS is 
not necessarily zero, P ZD of AMP-JS is plotted only for the M 
values in which P W zo < 0.1%. We see that with increasing the 
block length, the required number of measurements by AMP- 
JS is increased for a target Pzd- However, the performance of 
the proposed method does not depend on the block length. The 
comparison shows the superior performance of the proposed 
scheme, particularly for regular graphs with d v = 1. For very 
low number of measurements, the ltol graph outperforms the 
others. 

B. Noisy Measurements 

In Fig. [4j the probabilities of false alarm ( P FA ), detection 
(P D ), wrong zero detection ( P WZ d ) and zero detection (Pzd) 
are plotted versus the threshold c'. Both analytical and simula¬ 
tion results are presented. The sensing graphs are regular graphs 
with degrees d v = 1 and d M = 2. The number of blocks L is 
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Fig. 4: The probabilities Pd, Pfa, Pzd and Pwzd versus the 
threshold, c\ when the sensing graph is regular with degree distribution 
(1, 2) for L = 1000, a = 0.25 and SNR= 25 dB. The markers and 
the dashed lines correspond to the simulation and analytical results, 
respectively. 




Fig. 5: Pzd versus sparsity ratio for different noise levels with regular 
sensing graphs with the degree distribution of (1,2) for L = 1000 and 
M = 500. 


1000 and a = 0.25, and the curves are plotted for SNR= 25 
dB. As seen in this figure, the simulation and analytical results 
match closely. 

Fig-0 shows P ZD versus sparsity ratio for different noise 
levels when P W zd is limited below 2%. The curve for the 
noiseless scenario is also provided for reference. The sensing 
graphs are regular with degree distribution equal to (1,2), 
L = 1000 and M = 500. Fig. [5] shows that the performance 
improves with increasing SNR, and at each SNR, the P ZD drops 
rapidly when the sparsity ratio in increased beyond a certain 
threshold. Before such a threshold is reached, the performance 
is practically the same as that of the noiseless case, but, when 
this threshold is passed, the P ZD curve demonstrates a waterfall 
behavior down towards P zn = 0. This waterfall region (for 
a given SNR) corresponds to the values of a for which the 
threshold c' has to be decreased rapidly (with increasing a) to 
maintain P WZ d below 2%. As a consequence, P Z d is decreased 
sharply. 


VII. Conclusion 

A novel zero-block detection scheme in the context of wide¬ 
band spectrum sensing was proposed. For this, an analog to in¬ 
formation converter (AIC) with block sparse sensing matrix was 


designed. Analytical and simulation results for both noiseless 
and noisy scenarios were presented. The results demonstrated 
the effectiveness of the proposed scheme in reliable detection 
of spectrum holes with minimal complexity even in scenarios 
where the number of measurements were relatively small. 
In particular, it was shown that, at the presence of noise, 
performance similar to the noiseless case can be achieved at 
higher SNR values and smaller sparsity ratios. Probably one of 
the most important contributions of this work was to introduce 
a new compressed sensing paradigm in which one is interested 
in reliable detection of (some of the) zero blocks rather than 
the recovery of the whole block sparse signal. Future work 
includes the optimization of irregular degree distributions for 
sensing matrix (graph) to maximize P D for a given P WZ n, 
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